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ABSTRACT 

We use hydrodynamical A^-body simulations to study the kinetic Sunyaev-Zel'dovich 
effect. We construct sets of maps, one square degree in size, in three different cosmo- 
logical models. We confirm earlier calculations that on the scales studied the kinetic 
effect is much smaller than the thermal (except close to the thermal null point), with 
an rms dispersion smaller by about a factor five in the Rayleigh- Jeans region. We 
study the redshift dependence of the rms distortion and the pixel distribution at the 
present epoch. We compute the angular power spectra of the maps, including their 
redshift dependence, and compare them with the thermal Sunyaev-Zel'dovich effect 
and with the expected cosmic microwave background anisotropy spectrum as well as 
with determinations by other authors. We correlate the kinetic effect with the thermal 
effect both pixel-by-pixel and for identified thermal sources in the maps to assess the 
extent to which the kinetic effect is enhanced in locations of strong thermal signal. 

Key words: galaxies: clusters, cosmic microwave background 



1 INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 
1972, 1980; for reviews see Rephaeli 1995 and Birkinshaw 
1999) is the change in energy experienced by cosmic mi- 
crowave background photons when they scatter from inter- 
vening gas, especially that in galaxy clusters. The dominant 
version from clusters is the thermal SZ effect, the gain in 
energy acquired from the thermal motion of the gas which 
is commonly at a temperature of tens of millions of degrees 
in clusters. The kinetic SZ effect is the Doppler shift arising 
from the bulk motion of the gas. 

The thermal effect has been quite well studied theoret- 
ically, and recently has become a burgeoning area of obser- 
vational activity with the construction of two-dimensional 
maps of clusters becoming commonplace (Jones et al. 1993; 
Myers et al. 1997; Carlstrom et al. 2000). By contrast, less 
attention has been given to the kinetic effect. It presents 
a considerable observational challenge, because it is pre- 
dicted to be much smaller than the thermal effect on the 
angular scales explored so far, and also because unlike the 
thermal effect it possesses no characteristic spectral signa- 
ture allowing it to be distinguished from primary cosmic 
microwave background anisotropies. It is also more diffi- 
cult to make semi-analytic calculations. The thermal ef- 
fect can be obtained fairly directly via the Press-Schechter 



(1974) approach; various aspects including number counts, 
the global CMB distortion, and its impact as secondary 
anisotropy for CMB measurements at small scales have been 
extensively studied (Cole & Kaiser 1988; Colafrancesco et 
al. 1994, 1997; Bartlett & Silk 1994; Barbosa et al. 1996; 
Eke, Cole & Frenk 1996; Aghanim et al. 1997; Komatsu & 
Kitayama 1999; Atrio-Barandela & Miicket 1999; Molnar & 
Birkinshaw 2000; Hernandez-Monteagudo, Atrio-Barandela 
& Miicket 2000). Predictions for the kinetic effect require 
the simultaneous estimation of both mass and peculiar ve- 
locity, and the signal is much less dominated by the gas 
which happens to be in massive halos. A detailed analytical 
model has recently been constructed by Valageas, Balbi & 
Silk (2000) [see also Benson et al. 2000], which computes the 
effect including inhomogeneities in reionization; this super- 
sedes earlier modelling by Aghanim et al. (1997) where the 
peculiar velocities were drawn from a gaussian distribution 
of fixed width. 

We recently used large-scale hydrodynamical TV-body 
simulations to make simulated maps of the thermal effect 
and analyzed their properties (da Silva et al. 2000; for sim- 
ilar recent work see Refregier et al. 2000a; Seljak, Burwell 
& Pen 2000; Springel, White & Hernquist 2000). In this 
paper, we make a detailed analysis of maps of the kinetic 
effect, made using the same technique. Recently Springel 
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et al. (2000) used similar simulations to study the angu- 
lar power spectrum of the kinetic effect in the ACDM cos- 
mology, and investigated the influence of non-gravitational 
heating upon it. In this paper we study three different cos- 
mologies, and investigate a wide range of properties of the 
resulting maps including the relation between thermal and 
kinetic distortion. 



2 THE KINETIC SZ EFFECT 

2.1 Linear and non-linecir effects 

Traditionally, the phrase "kinetic SZ effect" is associated 
with the scattering of photons from galaxy clusters or other 
large gravitationally-bound objects. However, in simulations 
there is no clear distinction between this and the scattering 
from any other moving gas, and when we construct maps 
they will be based on the bulk motions of all the mate- 
rial in the simulation box. However, one should recognize 
that part of this signal is computable in linear perturba- 
tion theory, and is for example computed in a run of the 
CMBFAST program (SeljaJj & Zaldarriaga 1996) for a cos- 
mology with reionization, where it is called the low-redshift 
Dopplcr effect. These different regimes, and their relation 
to the Ostriker-Vishniac effect (Ostriker & Vishniac 1986), 
have been discussed in detailed by Hu (2000) and by Gnedin 
& Jaffe (2000), and observational aspects are described by 
Cooray & Hu (2000). 

Concerning the contribution from non-linear evolution, 
one should recognize that while there are non-linear con- 
tributions to the velocities, of particular importance is the 
non-linear evolution of the density field. This concentrates 
the matter into a virialized region which gives a high op- 
tical depth probe of the bulk velocity at that point. The 
non-linear effect is more important at low redshifts. 



2.2 The reionization epoch 

Photons only scatter at epochs where the gas is ionized. 
Reionization is thought to be induced by radiation from 
the first generation of massive stars in dwarf gala^xies (see 
Haiman & Knox 1999 for a review) , which unfortunately are 
well below our numerical resolution. Bruscoli et al. (1999) 
and Gnedin & Jaffe (2000) have modelled reionization using 
very small-scale high-resolution simulations, probing a range 
of angular scales complementary to our own, and compute 
the corresponding microwave anisotropics. We will not con- 
sider inhomogoneous reionization in this paper (sec Aghanim 
et al. 1996; Grusinov & Hu 1998; Knox, Scoccimarro & Do- 
delson 1998; Benson et al. 2000; Valageas et al. 2000), and 
we will treat the epoch of reionization as an independent pa- 
rameter. According to Gnedin & Jaffe (2000), anisotropics 
from inhomogeneous reionization will be subdominant on 
the scales we consider, especially if massive stars arc indeed 
responsible for reionization. 

For the thermal SZ effect the epoch of reionization is ir- 
relevant, as significant signals only come from low redshifts 
and from regions hot enough to be fully ionized. This is no 
longer true for the kinetic effect, where the increased den- 
sity of material at high redshift (per observed solid angle) 
maintains a significant signal even though the velocities are 



smaller. Indeed, we will see that, unlike the thermal SZ sig- 
nal, the kinetic SZ signal is not convergent and results are 
mildly dependent on the assumed redshift of reionization 
(though this is certainly not the largest uncertainty in our 
calculations). For case of comparison between cosmologies, 
we will assume that reionization took place at redshift 10, 
which satisfies all existing observations (see e.g. Griffiths 
et al. 1999; Haiman & Knox 1999), independent of the cos- 
mology even though detailed calculations would favour later 
reionization in the high-density case. 



2.3 The equations 

The thermal SZ effect in a given direction is computed as a 

lino integral, which gives the Compton j/-parameter as (Sun- 
yaev & Zol'dovich 1972, 1980) 



meC^ 



Tb Tie dl . 



(1) 



In this expression To and rio are the temperature and num- 
ber density of the electrons, the Boltzmann constant, ctt 
the Thomson cross-section, c the speed of light and me the 
electron rest mass. The frequency dependence of the effect is 
well documented (e.g. Birkinshaw 1999); ignoring relativis- 
tic corrections, it gives a temperature fluctuation 
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where x = Iw/kjiT is the dimcnsionless frequency and 
T the radiation temperature. In the long-wavelength limit 
X <C 1 (the Rayleigh- Jeans portion of the spectrum) we have 
ATth/T ~ — 2j/, and the thermal effect has a null point at 
X ~ 3.83, corresponding to ~ 217 GHz. 

Whenever the gas has a peculiar velocity in the CMB 
radiation frame, there is a kinematic Doppler effect in ad- 
dition to heating from the random thermal motions of the 
electrons. The change in the spectral intensity due to the 
kinetic effect is (Sunyacv & Zel'dovich 1980) 
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where 7o = 2kl,T^/h^c^, Vr is the line-of-sight velocity of the 
gas (positive/negative when receding/approaching) , h{x) is 
given by 



h{x) 



(e- - 1)'' ' 
and r is the optical depth 
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The temperature fluctuation due to the kinetic effect is 

■ (7t Tie dl , (6) 



which, in contrast to the thermal effect, is a frequency- 
independent quantity. As such, the kinetic effect cannot be 
distinguished from primary microwave anisotropics on the 
basis of its spectral dependence. 
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3 THE MAPS 

For full details as to the construction of maps, see da Silva 
et al. (2000). The starting ingredient is a set of three hydro- 
dynamical simulations, one for each of three cosmologies. 

• ACDM: a low-density model with a flat spatial geome- 
try and r^o = 0.35 and f^A = 0.65. 

• rCDM: a critical-density model with S7o = 1 and 
= 0. 

• OCDM: a low-density open model with Q,o — 0.35 and 
Oa = 0. 

The simulations were carried out using the public domain 
Hydra code (Adaptive P^M-SPH: Couchman, Thomas & 
Pearce 1995). In each case the box-size was 150 Mpc, 
and equal numbers of dark matter and gas particles were 
used. 

To construct a map we stack simulation output boxes, 
with random translations and reorientations, to reach high 
redshift. Outputs from the simulations were made at dif- 
ferent epochs so that they can be stacked smoothly into 
an evolving sequence. Typically around 35 boxes need to 
be stacked to reach to redshift 10; however at the redshifts 
from which the bulk of the signal arises only a small frac- 
tion of the volume of the box contributes to the maps and 
so we are not limited by cosmic variance in having only a 
single hydrodynamical simulation. For each cosmology we 
make 30 different maps, corresponding to different random 
orientations of the stacked simulations. 

To compute the thermal and kinetic effects we follow 
the approach described in da Silva et al. (2000) . In the sim- 
ulations each gas particle occupies a volume with radius pro- 
portional to its SPH smoothing radius, hi, inside which the 
mass profile is taken as mgasVK(r — ri,/ii), where n is the 
position of the particle's centre, mgas is the mass of the 
gas particle and W is the normalized spherically-symmetric 
smoothing kernel. 

The quantities we are interested in mapping are y and 
ATk/T. To do this we first convert the line-of-sight inte- 
grations in equations (|l]) and (^) into volume integrals over 
space. The resulting expressions are then discretized into cu- 
bic cells of volume V (voxels), which are stacked along the 
line of sight to produce pixels of area A. The total tempera- 
ture distortion induced by the kinetic effect in each pixel of 
the map can then be evaluated as 

An ax V 0.88 mga, 

— = --A^;^Z.Z.^-W/(|r.-n|,fe), (7) 

i CK 

and likewise for the thermal SZ y-parameter 

i ct 

where the i index runs over all particles which contribute to 
the pixel columns, and the a index is over the line-of-sight 
cubes. For each particle the sum of the smoothing kernel over 
voxels is normalized to 1/V. The quantities y and ATk/r 
are the average values in each pixel. In these expressions, 
Vi^i and Ti are the velocity and the temperature of the gas 
particles, rup is the proton mass, and the factor 0.88 gives 
the number of electrons per baryon, assuming a 24 per cent 
helium fraction and complete ionization in the regions of in- 



terest. We construct maps of the individual simulation boxes 
in our stacking, which are then added to give the final maps. 



4 RESULTS 

Typical example maps in each cosmology are shown in Fig- 
ure 0, with the same colour scale in each.Q The original 
maps have been smoothed using a gaussian with a full-width 
half-maximum (FWHM) of 1'. Unlike the thermal maps, the 
kinetic maps show positive and negative temperature fluc- 
tuations. The low-density maps (top and bottom panels) 
present stronger fluctuations than the critical density ones 
(centre panels). The brightest objects in the thermal maps 
do not necessarily correspond to the largest kinetic distor- 
tions (and vice- versa). However, bright thermal sources often 
reveal neighboring kinetic peaks of opposite signs (Diaferio, 
Sunyaev & Nusser 2000), resulting from mergings of sub- 
structures or rotational flows within the sources. This can 
be seen, for example, corresponding to the two brightest ob- 
jects of the thermal ACDM map and in the bright source on 
the lower right corner of the OCDM map. These objects are 
at redshifts z=1.7, 1.8, and 1.9 respectively. Fainter super- 
clusters at lower redshifts also produce similar signatures 
(see animations of individual maps on our WWW page), 
but their contribution to the final map tends to be erased 
by the box stacking process. Although the kinetic effect has 
no spectral signature, the search for these characteristic pos- 
itive/negative features in the SZ flux of bright sources and 
nearby super-clusters may prove to be an additional aid in 
isolating the kinetic effect in future SZ measurements and 
in identifying super-clusters at cosmological distances (see 
Diaferio et al. 2000). 

4.1 The redshift dependence 

In order to test how well the calculation of the effects con- 
verges with redshift. Figure ^ shows the rms distortion, 
(Jk — \J {{An/TY) , in the 1 arcminute resolution maps 
constructed out to redshift z, as a function of that redshift. 
Unlike the dispersion of the thermal effect, shown in Fig- 
ure ^, the integrated kinetic signal is not strongly conver- 
gent to high redshift, with the smaller velocities being com- 
pensated by the increased amount of material contributing. 
In the real Universe the signal is cut off by the epoch of 
reionization, and as can be seen our results are modestly 
dependent on that choice (though other uncertainties such 
as the normalization of the matter power spectrum will be 
more signiflcant). As advertised in Section 2, we choose to 
construct maps with data out to redshift 10 in each cosmol- 
ogy- 

4.2 The dispersion 

In Table 1 we list the rms fluctuations in both the thermal 
and kinetic effects, for a range of different smoothings to rep- 
resent possible instrumental resolutions. They are given as 

* A more extensive selection of colour maps, along witli anima- 
tions showing the contribution from each redshift, can be found 
at star-www. cpes . susx . ac.uk/'^andrewl/ sz/ sz .html 
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Figure 1. Pairs of simulated SZ maps of size one square degree. The left panels show AT^ij/T in the Rayleigh— Jeans part of the 
spectrum, while the right panels show the corresponding kinetic maps of ATjj/r from the same map realization. In each case, the original 
map was convolved with a gaussian beam profile of FWHM = 1'. From top to bottom, the three sets are ACDM, rCDM and OCDM. 
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Figure 2. The rms kinetic distortion in the 1' resolution maps 
made from individual simulation boxes at a redshift z (upper 
panel), and the rms of the accumulation of these maps up to 
redshift z (lower panel) . Each curve is an average over the 30 map 
realizations. Note that the lower panel is not simply the integral 
of the upper, because when kinetic maps are added together they 
'interfere'. 



Table 1. The rms fluctuations of the thermal (c!t\i) and kinetic 
((T]^) SZ effects as a function of beam resolution in the Rayleigh- 
Jeans limit. Values are in multiples of 10"'^. 
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temperature fluctuations in the long-wavelength part of the 
spectrum. Because the thermal effect is purely additive, its 
fluctuations, which is all that most instruments are capable 
of measuring, do not dominate over the kinetic fluctuations 
by as large a factor as one would surmise from looking at 
the total magnitude of the effects. 



4.3 The pixel distribution 

Figure ^ shows a histogram of the pixel distribution in the 
kinetic maps for the ACDM cosmology. The different curves 
represent different smoothings ranging from 0.5' to 10'. The 
curves show the distribution of ATk/T values which would 



Figure 3. As Figure ji] but for the thermal effect. Strong con- 
vergence to high redshift is seen. Note that this plot, showing 
the dispersion, differs from a similar plot in da Silva et al. (2000) 
which showed the convergence of the mean y-distortion with red- 
shift. 



be seen if an instrument with the corresponding resolution 
points randomly at the sky. At low resolutions (5 arcminutes 
or worse) the curves are close to gaussian distributed for all 
intensities, whereas at higher resolution marked nongaussian 
tails are present from non-linear structures. These nonlinear 
tails represent the best hope for observing the effect directly. 
As expected the pixel distribution mean values are very close 
to zero. 

These pixel distributions are quite discouraging for at- 
tempts to detect the kinetic effect from individual sources 
with the Planck satellite, whose best resolution is about five 
arcminutes. We see that the pixel distribution at such reso- 
lution is well approximated by a gaussian even for the rarest 
events we see, with no non-gaussian tail of bright rare events. 
The rms dispersion is very small compared to the expected 
level of primary anisotropics on those scales, and so the ki- 
netic effect will not be detectable directly without some ex- 
ternal information as to where to look to see it. Later in 
this paper, we shall consider the use of the thermal effect 
as an indicator of where to look for the kinetic effect. This 
can allow a statistical measure of the cluster velocity field 
(Aghanim et al. 1997; see also Ifaehnelt & Tegmark 1996; 
Ifobson et al. 1998). Indeed, the MAP and Planck missions 
could allow the overall detection of the cosmic bulk flow 
when coupled with X-ray observations (Kashlinsky & Atrio- 
Barandela 2000). 
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Figure 4. A histogram of the AT]^/T values in the ACDM maps 
with different smoothings. The dashed lines show gaussian fits to 
the distributions. 



4.4 The angular power spectrum 

We calculate the power spectrum using the flat-sky approxi- 
mation in the usual way; see e.g. Peacock (1999). In Figure ^ 
we show the thermal and kinetic SZ angular power spectra 
obtained from simulations for the three cosmologies, with 
the Ci band averaged into logarithmic bins in the Fourier 
space. 

The power spectra shown are averages over the thirty 
maps we made for each cosmology. To estimate the statisti- 
cal error we used bootstrap resampling of the thirty maps, 
which gave a range of estimates indicated at 68 per cent 
confidence by the dotted lines (note that the map-to-map 
variations are of course much bigger than the variation in 
the bootstrap resamples of thirty maps each). For the two 
low-density cosmologies we see that the statistical errors 
are small, while for the critical-density thermal case even 
thirty maps leaves a significant residual error, because the 
maps tend to be dominated by the presense or absence of 
single bright features. Unfortunately computer time limita- 
tions prevent us increasing the number of maps used to make 
the estimate, and from running more sets of hydrodynami- 
cal simulations to remove possible correlations from using a 
single set. 

From Figure ^ we see that the low-density models 
present similar shapes and amplitudes for the thermal and 
kinetic angular power spectra. In the critical-density case 
both effects produce less power at all scales, which can be 
explained by there being fewer objects along the line of 
sight for this model. On large angular scales, the thermal Ci 
curves typically exceed the kinetic ones by a factor of about 
25 for all cosmologies. On smaller scales {£ ^ 2000) the gap 
between the thermal and kinetic power spectrum increases 
slowly in the low-density models, whereas in the rCDM 
model it decreases rapidly as the thermal Ce curve falls 
with £. This can be understood if we take into account that 
although the simulations are normalized to give the same 
abundance of objects at redshift zero, in the low-density 
models there are SZ sources visible to much higher redshifts. 
As we will see below, on small angular scales the thermal SZ 
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Figure 5. In each panel, we see the primary CMB, thermal SZ 
and kinetic SZ angular power spectra obtained for each cosmol- 
ogy. From top to bottom we have ACDM, rCDM and OCDM, 
respectively. The thermal spectrum is plotted for the Rayleigh— 
Jeans part of the frequency spectrum, and can be scaled to other 
frequencies. The dotted lines indicate 1-sigma errors, obtained by 
bootstrap resampling the thirty map realizations. 



spectrum receives strong contributions from sources above 
redshift one. The lack of distant sources in the rCDM model 
is then reflected in the thermal SZ power spectrum by a 
strong reduction of power on these scales. 

In the ACDM cosmology the amplitude of the thermal 
SZ power spectrum becomes comparable to the power spec- 
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Figure 6. The redshift dependence of the power spectrum of the 
kinetic effect is shown in the upper panel, with the lower panel 
showing the thermal effect for comparison, fn each case, the spec- 
trum is shown for the accumulated signal from redshifts greater 
than a given value. These plots refer to the ACDM simulation. 

trum of the primary CMB around £ = 2000. This prediction 
agrees with earUer resuhs from numerical simulations (Re- 
fregier et al. 2000a; Seljak et al. 2000; Springel et al. 2000) 
and confirms that the SZ power spectrum should not be 
an important source of contamination for MAP (its ampli- 
tude is well below the projected noise for this mission, see 
Refregier et al. 2000a, b), whereas it should be detectable 
by Planck and other experiments probing smaller angular 
scales 

Figure ^ shows the dependence of the thermal and ki- 
netic angular power spectra with redshift. We see that the 
kinetic effect has had a much weaker dependence on redshift 
(for £ ~ 2000 about 15 per cent of the power has origin above 
z > 7.5), and indeed, while sub-dominant now, was domi- 
nant at redshifts earlier than about three. On large angular 
scales, the thermal power spectrum is clearly dominated by 
sources at low redshifts. More than two-thirds of the power 
on these scales is generated at redshifts z < 1, whereas on 
small scales the majority of the signal is produced above 
redshift 1. On intermediate scales there is still a reasonable 
fraction of power with origin above z — 1; for instance at 
£ = 2000 that fraction is about 50 per cent. 

In general our results show broad agreement with the 
SZ power spectra obtained from numerical simulations by 



Figure 7. A scatter plot of the kinetic effect against the thermal 
effect, pixel-by-pixel, for a single ACDM map. 

Refregier et al. (2000a), Seljak et al. (2000) and Springel et 
al. (2000). In the first two cases the authors have also studied 
three CDM models, but made no predictions for the kinetic 
SZ power spectrum. In Springel et al. (2000) the angular 
power spectrum is computed for both effects, but results are 
only reported for the ACDM cosmology. After taking into 
account the difference between simulations on as and Jls/i 
we find that the thermal spectra agree quite well on large 
angular scales. In particular, our thermal ACDM spectrum 
agrees remarkably well with Springel et al. (2000) for a wide 
range of £ (the difference is within 20 per cent for £ ^1200). 
This is an encouraging result as all simulations should agree 
on these scales (as shown in Refregier et al. 2000a, the finite 
size of the simulation boxes has very little influence on the 
SZ power spectrum). 

On small angular scales different simulations give rather 
different results for the thermal effect. Generally speaking 
our simulations predict less power than Springel et al. (2000) 
and more power than Refregier et al. (2000a) and Seljak et 
al. (2000). However, in the critical-density case our ther- 
mal curve peaks before and falls faster than the lower limit 
for the spectrum reported by Refregier et al. (2000a). The 
fact that Springel et al. (2000) predict more power on these 
scales is most likely reflecting the higher resolution of their 
simulations. 

The only available results for comparison of the kinetic 
spectra are those of Springel et al. (2000) for the ACDM cos- 
mology. The agreement between simulations is not as good 
as in the thermal case, as the shapes of the kinetic curves 
differ on large scales. Our simulations predict about ~ 3.5 
times less power at ^ ~ 400. However simulations show good 
agreement in the range 1000 ^ £ ^ 2000. Different assump- 
tions concerning the ionization history of the gas may be the 
main source of differences in the spectra at very high £. 

4.5 Kinetic effect versus thermal effect 

Given the smallness of the kinetic effect, it is advantageous 
to have an idea where on the sky to look in order to see it. 
Such a guide is given by the thermal effect, which indicates 
locations where gravitational collapse has concentrated scat- 
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Figure 8. For a given value of thermal distortion, this plot shows 
the dispersion of the values for the kinetic distortion found at the 
same spatial locations. For the plotted values, the distribution is 
quite well approximated by a gaussian, whereas for larger values 
it is dominated by nongaussian effects. 



tering electrons. Figure |7| is a pixel-by-pixel scatter plot of 
the kinetic effect versus the thermal effect, found in a single 
map of tlie ACDM simulation. The kinetic signal is roughly 
symmetric around zero and its dispersion increases for the 
pixels with higher thermal distortions. The 'tendrils' are the 
pixels associated with particular large bright features in the 
maps. 

An interesting question is what values of kinetic distor- 
tion are expected at locations with a given thermal distor- 
tion, and this is shown in Figure ^. We binned the pixels 
of all thirty ACDM maps according to their thermal values 
and studied the distribution of kinetic distortions. As long 
as the distortion is not too large, these distributions are well 
approximated by gaussians of width (Jk as shown in the fig- 
ure, the brighter thermal sources corresponding typically to 
larger kinetic distortions. This enhancement at locations of 
larger thermal signal is what allows a statistical detection 
of the magnitude of the kinetic effect (Haehnelt & Tegmark 
1996; Aghanim et al. 1997) and may allow detections in 
some rare bright objects. For a discussion of identifying the 
kinetic effect through the component separation process, see 
Hobson et al. (1998). 

A similar exercise can be carried out for identified 
sources rather than pixels. To mimic typical observational 
procedures, we first subtract the mean signal from the ther- 
mal maps. We then identify source locations in the ther- 
mal maps using the SExtractor package (Berlin & Arnouts 
1996), and determine the total thermal flux of the sources 
(at a given frequency channel) by adding up the flux in all 
contiguous pixels until the pixel intensity falls below half 
the maximum. In order to reduce source confusion effects 
we only consider bright resolved sources which satisfy the 
selection criterion of having a half maximum intensity above 
twice the rms distortion in the maps. We then compute the 
corresponding flux of the kinetic effect by summing the same 
set of pixels in the kinetic maps. Figure H shows the resulting 
scatter plot, at frequency 143 GHz, made from all thirty of 
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Figure 9. As Figure ^, but now showing the total flux of sources 
as described in text. 



the one arcminute resolution ACDM maps (unlike the pixel 
scatter plots which are from a single map). 

As expected the largest kinetic distortions correspond 
to sources with stronger thermal fluxes. However, at this 
frequency the typical flux of the kinetic effect is around one 
order of magnitude below the thermal, though for specific 
clusters the ratio can be as high as one third. 



5 CONCLUSIONS 

We have used hydrodynamical simulations to study a num- 
ber of aspects of the kinetic SZ effect, including its depen- 
dence on cosmological parameters. We have studied the red- 
shift dependence, pixel histograms, angular power spectra, 
and the correlation between the kinetic and thermal effects. 
We have confirmed that the kinetic effect has a dispersion 
more than a factor five below the thermal (in the Rayleigh- 
Jeans region) , leading to an angular power spectrum a factor 
of typically twenty-five lower in the low-density cosmologies. 
For critical density we have found a smaller difference, but 
with much greater statistical uncertainty. For the thermal 
effect the angular power spectrum is mostly generated at 
redshifts below one, while for the kinetic effect a significant 
amount of the power has origin above redshift two. The cor- 
relation of the kinetic effect with the thermal confirms and 
quantifies the expected enhancement of the kinetic effect in 
regions with a strong thermal signal. 

Our simulations are of ideal size to study the SZ ef- 
fect on scales between one and several arcminutes, which is 
currently a resolution attracting great interest. Our results 
complement perfectly the small-scale simulation work of Br- 
uscoli et al. (1999) and Gnedin & Jaffe (2000), and recent 
semi-analytic work including that of Benson et al. (2000) 
and Valageas et al. (2000). The theory of these secondary 
anisotropies is now becoming highly developed; the key chal- 
lenges for the kinetic effect lie very much on the observa- 
tional side. 
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